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Reheating after inflation can occur through inflaton decay or efficient parametric resonant produc- 
tion of particles from the oscillation of the inflaton. If the particles produced interact with scalars 
that were light during inflation, then significant super-horizon density perturbations are generated 
during this era. These perturbations can be highly non-Gaussian. 



Measurements of the cosmic microwave background ra- 
diation 0, U have clearly shown the presence of super- 
horizon primordial density fluctuations at roughly one 
part in 10^. Inflation provides a natural explanation for 
such density fluctuations, since vacuum fluctuations of 
the inflaton (or any other light scalar field) get pushed 
outside of the horizon and enter at a much later time 
as classical density perturbations W\. Recently, Dvali, 
Gruzinov and Zaldarriaga Jk, 5j and Kofman 6] (DGZK) 
have shown in a number of scenarios how the interac- 
tions of such additional light fields to, e.g. the inflaton, 
could also generate adiabatic density fluctuations, inde- 
pendent of those created by the inflaton dynamics. In 
this scenario the size of non-Gaussian perturbations can 
be much larger than what occurs in single-field inflation- 
ary models j^, ■ 

This is achieved by coupling a light scalar field to a 
heavier field that at some time subsequent to inflation 
dominates the energy of the Universe, such that the par- 
ticle properties of this heavier field are modified by the 
fluctuations of the light field. When the heavier particle 
decays, spatial fluctuations in either its mass or its de- 
cay width generate energy density perturbations in the 
radiation. This is because before reheating the universe 
is matter dominated, with the oscillating heavier particle 
dominating the total energy density, while after the de- 
cay the universe is radiation dominated. As the energy 
density in matter redshifts slower than energy density in 
radiation, regions of the universe where the decay occurs 
at a later time stay matter dominated longer and will be 
denser than regions where decay happens earlier. This 
gives density perturbations of order 
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is the time between the end of inflation [to) and reheating 
(^rh)- The evolution of density perturbations in this 
scenario has been studied in detail in 0, Q . 

In a similar way, modifications to the particle proper- 
tics of the particles produced during reheating can also 
introduce energy density perturbations. Density pertur- 
bations are created if the decay products interact with 
fields that were light during infiation. 



To see this, we need to discuss how the inflaton re- 
heats. Suppose reheating occurs through direct (Born) 
decay of the inflaton. Then a fluctuation in the mass of 
the decay product x modifles the inflaton decay width, 
because of the dependence of the available phase space 
on the masses of the final state particles. These lead to 
calculable density perturbations since the exact depen- 
dence of the width on the mass of the light particles can 
be computed in any given model Ifj for example, 

the infiaton decays via $ XX then the tree-level decay 
width is modified from phase space by an amount 
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If $ decays near threshold, then the resulting density 
perturbation dependence on (5to^ can be large and highly 
non-linear. 

We expect a (5m^ with a super-horizon spatial varia- 
tion to be generated if x interacts with a field a that was 
light during the infiationary era and through to the era 
of reheating. Note that even in the absence of direct cou- 
plings of the fields x and a, they are expected to interact 
indirectly through some intermediate states. Quantum 
corrections will typically generate a dependence of to^ 
on the super-horizon fiuctuations of a at some order in 
perturbation theory, as indicated by Fig. ^ In this pa- 
per we focus mainly on the effect that fiuctuations in 
the mass of the particles produced during reheating or 
preheating have on density perturbations. 

Besides reheating through direct Born decay, the infla- 
ton may instead reheat the universe through parametric 
resonance (preheating) |M lli ll^ lisl llj. Preheating 
can be very efficient and be completed very soon after 
inflation, within 0(10 — 100) oscillations of the inflaton 
fleld about its minimum. Whether this process of re- 
heating dominates over the Born decay into bosons or 
fermions depends on the parameters of the model ^ . 

In the scenario of DGZK, additional density pertur- 
bations can be created during preheating, by modifying 
the time it takes for parametric resonance to complete 



^ The growth of perturbations during the m atte r-dominated era of 
the oscillating inflaton has been studied in Il5l and, if parametric 
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FIG. 1; Quantum corrections may generate a dependence of 
m?, on super-horizon fluctuations in a. 



and for the universe to thermalize. The size of this time 
interval depends on the parameters of the model, and in 
particular on the mass of the produced particles, which 
we discuss below in a simple model. This is the main 
subject of this paper. Depending on how efficient pre- 
heating is, the size of the time interval can have a weak 
or strong sensitivity to the mass of the decay products. 

We use the canonical model of preheating and add a 
scalar a which we assume is light during inflation so that 
it acquires super- horizon perturbations 5a {x) ~ ^^inf 
during that era. For this to occur it is necessary that 
during inflation its mass satisfles mg. < -ffjnf. cr is as- 
sumed to interact more strongly with the x compared to 
The interactions we consider are 
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A Z2 symmetry x ~^ ~X been imposed for simplicity. 
Self-interactions cr^ and x^ are assumed to be irrelevant 
during the first stage of preheating defined below. We 
assume that at the end of inflation the fields x S'Hd ^'^^ 
near enough to the minimum of their potential so that 
we can neglect the motion of their zero modes. 

Inflation ends when t = to — and is followed by 

a matter-dominated era described by rapid oscillations 
of the inflaton about the minimum of its potential which 
we assume to be 
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For simplicity we assume that the inflationary potential is 
also described by this simple quadratic form, giving rise 
to chaotic inflation 0|. During inflation, _ffinf ~ 777,$. 
At the end of inflation $ = $0 — JTipi/S and thereafter 
decays as ^(t) ~ ■m^xjZm.^t. 

For large enough coupling these oscillations trig- 
ger parametric resonance, and the energy density in x 
increases exponentially yii]. If this process is efficient, 
the universe eventually is dominated by the x particles, 
which then thermalize the universe at some later time 
through its interactions with Standard Model or Grand 
Unifled Model particles. 

The perturbations in the inflaton give rise to adiabatic 
density perturbations, whose size depend on the form of 
the inflaton potential. In this letter we concentrate on the 
density perturbations generated from the fluctuations in 



the cr scalar field. In de Sitter space [Tsjl 
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where in Eq. is the number of e-foldings during 
infiation. In Eq. |(7J) the comoving coordinates x, y are 
well seperated and we neglect the logarithmic dependence 
on \x — y\. 

The X field does not acquire super-horizon perturba- 
tions because its effective mass 
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during inflation is larger than the Hubble parameter for 
parameter values which allow for efficient parametric res- 
onance. Henceforth we absorb the XH^^^N/ (Air^) into 
TO^ . Treating a as an external field, its fluctuations can 
be absorbed into fiuctuations in the mass of the field %, 

Sml = 2fiS(7 + X5a^. (9) 

where we have used 

(10) 



Sa — a — (a) , Sa'^ = cr^ — (cr^) 



and we will impose (cr) = 0. The size of the fluctuations 
Srn^ is determined by the two-point function 

{Sml{x)6ml{y))=4f^^<j{x)a{y))+2X^{<jix)a{y))' 

Using Eqs. © and Q we find the fluctuations for widely 
separated comoving coordinates x and y to be of order 
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While the field cr(a;) is Gaussian, the fiuctuation Srn^ 
is only Gaussian for A = 0. For XH^^^ » fJ-Hinf, Snv^ is 
highly non-Gaussian. For example, consider in this limit 
the three-point function for equally separated comoving 
coordinates. One finds for the analog of skewness 



{S'm'^{x)Sm'^{y)S'm'^{z)) 
(Jm^ (x)(5m2 (y))3/2 

_ 8{<j{x)aiy)){a{y)<7{z)){a{z)a{x)) 
23/2(a(a;)cr(y))3 



(12) 
= 2V2. 



For the remainder of this paper we set /i = which cor- 
responds to imposing a a —a symmetry. We make 
this decision to simplify the analysis of the backreaction 
of X on cr discussed below. Then 
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These fluctuations in Sra^ are non-Gaussian and always 
positive. 
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This situation would be excluded if this were the only 
source of density perturbations. A more interesting sce- 
nario in this situation would be if the dominant source 
of perturbations came from the inflaton potential. Then 
the perturbations generated during preheating providing 
a sub-dominant, non-Gaussian contribution. Since here 
the source for the non-Gaussian perturbations is not the 
same as the source -the inflaton- providing the domi- 
nant Gaussian contribution, the current limits on non- 
Gaussianity | T9| do not apply, since those limits assume 
that the non-Gaussian and Gaussian perturbations are 
generated by the same field. 

We define preheating to last until significant particle 
production of x occurs and the energy densities in $ and 
X become equal. The duration of this stage depends on 
and coupling constant g, 



Fluctuations in and the coupling g give rise to density 
fluctuations from Eq. (Q. 

Fluctuations in g can be generated if it is replaced by 
an effective coupling 
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where M is some mass scale j2^]. The a dependence of 
(7eff generates non-Gaussian perturbations 5g = Sg"^ j g^ — 
Hf^^/M"^ . It also modifies the large time-dependent mass 
of X, an effect that is distinct from modifying m-^. Here 
too we have to worry about the backreaction of x on a. 

Next we describe our numerical method for determin- 
ing the energy density in x during preheating. Neglect- 
ing the backreaction of x on the inflaton, which only 
becomes significant at the end of the preheating stage 
when = pcj) 12], the equation of motion for the fields 
X = x(ao/a)'^/^ (a is the scale factor) and a are 

Xl' + [Ak + 2gcos(2(z - zo))] Xfc = (16) 



5a" + -Ja' + m^_„ff(5a = , (17) 

where derivatives are with respect to z = m^t and we 
have chosen zq = 1. We have defined 

q = ^ = qo^ , ao = a{to) 
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and the mass parameters are given by 



(19) 



The equation for Xk describes a time-dependent 
harmonic oscillator with frequency fl^ = 



m^iAk + 2qcos{2{z — Zq))]. In the limit of a static 
universe and constant Sa this equation reduces to the 
Mathieu equation. 

Efficient parametric resonance requires go S> 1 and 
fh^ < m^. Note that we included a term of order XH^^^N 
into the definition of m^, where N is the number of e- 
foldings during inflation. The bound ^ therefore 
implies N < A^^. For the values of A we consider, this 
is a very weak bound on the number of c-foldings during 
inflation. 

For a given value of k the energy density in x is 
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where Nk (t) is the number density for a mode with given 
wave number k. The number density can be calculated by 
numerically solving for the Bogolyubov coefficient, giving 
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with initial conditions Xfc(io) 



TO*Xfc(*o) 



rife/2. The field Xfe satisfies the same equation as Xfe 
and is related to it (see Appendix B of |3). The energy 
density is obtained by integrating Eq. (|2()|l to obtain 
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The exponentially large number density of x particles 
leads to a large backreaction on a that must be included 
to correctly determine the size of the effect we are de- 
scribing. The backreaction of x on can have two effects: 
first, it can lead to production of large numbers of a par- 
ticles, and second it gives rise to a large effective mass of 
the cr field. The first effect was analyzed by Felder and 
Kofman using a numerical lattice simulation of pre- 
heating and the subsequent thermalization of the x with 
the a fields. In their Figures 14 and 15 they show the 
number densities of $, x and cr. Their numerical results 
show that during preheating the number density in cr is 
much smaller than in cither x or $ and its effect on the 
evolution of either or is negligible. The second 
effect is more significant. Once mo-, off gets larger than 
H, the amplitude 6<7 will decrease rapidly {2^ . To sim- 
plify the analysis we will assume that the dependence of 
mcr,cff on iTicr CBii be neglected. To estimate the time at 
which the backreaction becomes important, we compare 
the effective mass m'^ ~ A(x^) to the Hubble parame- 
ter. The ratio that determines their relative importance 
can be expressed as 
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where we have used m^^off — g\^\ and ~ (7|<i>|np(, ~ 
.g2$2(^2^ El- For A - IQ-^ and H ~ 2m$/300 we find 
that this backreaction becomes important when p^/ ~ 
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3g^. For A ^ 10~^ the backreaction becomes important 
when Px/p<s> ~ 300 g^. In this letter we wiU not solve the 
full coupled set of differential equations, but rather deal 
with this backreaction by turning off Sm^ at the time Zc 



when 
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Although for different values of Srn^ the intercept time Zc 
is different, that difference is second order in Sm'^. It is 
then sufficient to use the Zc obtained by setting Sm^ — 0. 
If Eq. (|24|l intercepts R along a plateau corresponding to 
no particle production, then we make the conservative 
choice of cutting off the mass fluctuation at the location 
of the first intercept. 

In Fig. 121 we display a logarithmic plot of the ratio ^ 
R = {t) I p$ (t) together with Eq. for scenario 1 , as 
defined in Tabled In order to estimate the sensitivity of , 
we show in Fig.|31a magnification 



2rh 



l+m<i)T on (5m^, 



of the region where /^(zrh) = 1- We also show these plots 
for the three other scenarios defined in Table ^ (keeping 
m^/iripi ~ 10^*^ fixed). 
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0.1 
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6 X 10" 
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0.4 
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0.06 



TABLE I: Definition of the four choices parameter sets. Also 
shown are the numerical results for Zc and Km, as defined in 
Eqs and (HSJ- 

We are interested in the change in r generated by a 
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FIG. 2: Logarithmic plot of i? = Px/P* 
eters are g = A x 10~*, = 10~^ 
shown is Eq. ||^ with A = 10"'^. 
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FIG. 3: Logarithmic plot of the effect of the mass fiuctuation 
(5m^ on «RH. The chosen parameters are given in each figure. 
The solid line corresponds to 5m\ — 0, while the long and 



fluctuation drri^/r 
can be expressed as 



A. Since A is tiny, that change short dashed lines correspond to 5ra^/ra^ = 10 and 2 x 



10" 



respectively. 
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From Tabled we can see that the typical Km is 0(0.1 — 1). 

The reader may wonder why we are using larger val- 
ues for 6m^ that are not consistent with the A we choose. 
Since the perturbation 6t is linear in (5m^, the k^, ob- 
tained this way is unchanged if we were to use smaller 
values for Sni^ . The reason for this choice of Sm^ is that 
the plots are easier to read. We also repeat that Zc was 
determined with the correct A. 

One may also wonder why the presence of a (5m^ at 
early times has any effect at all, especially given that it 
only persists while R ^ lO^'^ — 10~^. Parametric res- 
onance is dramatic because of stimulated emission. So 
even if at earlier times the production of x particles is af- 
fected due to a non-zero Snr^, this will impact the much 
greater growth occurring at later times. A more detailed 
numerical simulation, including all the effects of backre- 
action and scattering, such as done in plj for preheating 
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FIG. 4: Logarithmic plot of the effect of the fluctuation in 
the coupling constant Sg^ on zna- The chosen parameters are 
given in each figure. The solid line corresponds to Sg^ — 0, 
while the long and short dashed lines correspond to dg'^/g'^ = 
10^'^ and 2 x 10~^, respectively. 



without a fluctuating a field, is needed to explore in de- 
tail the sensitivity of Sp/p to super-horizon fluctuations 
in a. 

Mathematically, the intuition expressed above may be 
expressed in the following way. The density in x is ap- 
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proximately given by 
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dt v{t) 



to 
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where A'^ is a prefactor that depends on the parameters 
of the model. Here v is a, characteristic exponent leading 
to exponential growth. We approximate its dependence 
on k as given by its value near A; ~ 0. The coefficient v 
also depends on m^, so 
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i5m^ as occurring from the exponential. This gives 
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This result has 0(1) agreement with our previous nu- 
merical computations. (Compare Zc/ zbm with Hm-) It 
illustrates that St/t is not suppressed by any very small 
numbers other than Srn^/rri^. 

We also explore the dependence of t on fluctuations in 
5eff |20l | . For non-zero particle number the interaction 
(|15|l introduces a backreaction of x on cr corresponding 
to an effective mass Ti^^g = g^$^(x^)/M^. As before, 



we cut off the fluctuation in g^s when m 
occurs when 



2 

CT.cff 



Px 



M2 
2 



to: 



3i?2. This 



(29) 



The fluctuation in ^off 
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Numerically we find that vi/vq ~ and is positive. 

A negative correlation is expected, since both the char- 
acteristic exponents of the Mathieu equation and the in- 
stability bands are the largest near the kinematic limit 
A = 2q, corresponding to = k = 0. Increasing to^ 
removes more instability bands from the available phase 
space. Using the approximate formula above, we can 
solve for the change in the reheat time due to a fluctu- 
ation (5to^, approximating all the dependence of 6t on For $„ 



gives rise to non-Gaussian density perturbations. 

In Fig. ^ we display the ratio R ior g = A x 10~^, 
TO$/TOpi = 10^^ and to^/to| = 0.1. We choose two val- 
ues of M that give 6g = 10^''^ and Sg = 10"''. According 
to H29|l . the fluctuation in gcff is cut off at Px/p<s> = 10~^ 
and 10"^, respectively, corresponding to Zc — 26 and 
Zc — 47. For both of these parameters we find that there 
is a large linear effect which we express as 



— — KgOg . 

10~^ we find Kq = 0.9 and for Sg 



10- 



(31) 



we 



find Kg = 1.4. As in the previous case, in obtaining our 
plots we used larger values of 6g to determine Kg. 

In conclusion, we have shown that during preheating, 
interactions of the "decay products" of the infiaton with 
other light scalar fields can give rise to super-horizon 
mass fluctuations in these decay products. These fluctu- 
ations will then give rise to density perturbations of the 
universe. Depending on the coupling of the decay prod- 
ucts of the inflaton to the light scalar fields, the dominant 
density perturbations generated from this effect will be 
either Gaussian or non-Gaussian. 

This work was supported by the Department of Energy, 
under the contract DE-FG03-92ER40701. 



[1] G. F. Smoot et ai, Astrophys. J. 396, LI (1992); 

C. L. Bennett et at, Basic Astrophys. J. 464, LI (1996). 
[2] C. L. Bennett et ai, Astrophys. J. Suppl. 148, 1 (2003); 

D. N. Spergel et al. Astrophys. J. Suppl. 148, 175 (2003). 
[3] For a review see e.g. E. W. Kolb and M. S. Turner, The 

Early Universe, Perseus Publishing, (1990), Cambridge, 
Massachusetts. 

[4] G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 

69, 023505 (2004). 
[5] G. Dvali, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 

69, 083505 (2004). 
[6] L. Kofman, plenary talk Cosmo03, Chicago, September 

2002; astro-ph/0303614 
[7] M. Zaldarriaga, Phys. Rev. D 69, 043508 (2004). 
[8] S. Tsujikawa, Phys. Rev. D 68, 083510 (2003); 

A. Mazumdar and M. Postma, Phys. Lett. B 573, 5 

(2003), Erratum-ibid. B 585, 295 (2004); F. Vernizzi, 

Phys. Rev. D 69, 083526 (2004). 
[9] M. Postma, JCAP 0403, 006 (2004). 



[10] J. H. Traschen and R. H. Brandenberger, Phys. Rev. 
D 42, 2491 (1990); L. Kofman, A. D. Linde and 
A. A. Starobinsky, Phys. Rev. Lett. 73, 3195 (1994). 

[11] Y. Shtanov, J. H. Traschen and R. H. Brandenberger, 
Phys. Rev. D 51, 5438 (1995). 

[12] L. Kofman, A. Linde and A. Starobinsky, Phys. Rev. D 

56, 3258 (1997). 

[13] I. Zlatev, G. Huey and P. J. Steinhardt, Phys. Rev. D 

57, 2152 (1998). 

[14] R. Micha and I. Tkachev, Phys. Rev. D 70, 043538 
(2004). 

[15] B. Basset, Phys. Rev. D 56, 3429 (1997); H. Kodama 
and T. Hamazaki, Prog. Theor. Phys. 96, 949 (1996); 
T. Hamazaki and H. Kodama, Prog. Theor. Phys. 96, 
1123 (1996); Y. Nambu and A. Taruya, Prog. Theor. 
Phys. 97, 83 (1997); F. FineUi and R. Brandenberger, 
Phys. Rev. Lett. 82, 1362 (1999). 

[16] S. Khlebnikov and I. Tkachev, Phys. Lett. B 390, 80 
(1997); S. Khlebnikov and I. Tkachev, Phys. Rev. Lett. 



6 



79, 1607 (1997); S. Khlebnikov and I. Tkachev, Phys. 
Rev. D 56, 653 (1997); B. Bassett, D. Kaiser and 
R. Maartens, Phys. Lett. B 455, 84 (1999); B. A. Bas- 
sett, F. Tamburini, D. I. Kaiser and R. Maartens, Nucl. 
Phys. B 561, 188 (1999); F. FineUi and R. Branden- 
berger, Phys. Rev. D 62, 083502 (2000); F. FineUi and 
S. Khlebnikov, Phys. Lett. B 504, 309 (2001); F. FineUi 
and S. Khlebnikov, Phys. Rev. D 65, 043505 (2002); K. 
Enqvist, A. Jokinen, A. Mazumdar, T. Multamaki, and 
A. Vaihkonen, astro-ph/0411394 
[17] A. D. Linde, Phys. Lett. B 108, 389 (1982). 



[18] A. D. Linde, Particle Physics and Inflationary Cosmol- 
ogy, Harwood Academic PubUshers, (1990), New York, 
New York. 

[19] E. Komatsu et ai, Astrophys. J. Suppl. 148, 119 (2003). 
[20] F. Bernardeau, L. Kofman and J. P. Uzan, Phys. Rev. D 

70, 083004 (2004). 
[21] G. N. Felder and L. Kofman, Phys. Rev. D 63, 103503 

(2001). 

[22] K. Enqvist, A. Mazumdar and M. Postma, Phys. Rev. D 
67, 121303 (2003). 



